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Abstract: We predict level degeneracy of the rotational type in diatomic molecules described by means of a 
cotangent-hindered rigid rotator. The problem is shown to be exactly solvable in terms of non-classical 
Romanovski polynomials. The energies of such a system are linear combinations of t(t + 1) and 
l/[t(t+ 1) + 1/4] terms with the non-negative integer principal quantum number t = n+ \m\ being the sum of 
the order n of the polynomials and the absolute value, |m|, of the square root of the separation constant between 
the polar and azimuthal angular motions. The latter obeys with respect to t same branching rule, |m| = 0, 1, t, 
as does the magnetic quantum number with respect to the angular momentum, I, and in this fashion the t quantum 
number presents itself formally indistinguishable from I, In effect, the spectrum of the hindered rotator has the 
same (2t + l)-fold level multiplicity as the unperturbed one. For low t values wave functions and excitation 
energies of the perturbed rotator differ from the ordinary spherical harmonics, and the 1(1 + 1) law, respectively, 
while approaching them asymptotically with the increase of t. In this fashion the breaking of the rotational 
symmetry at the level of the representation functions is opaqued by the level degeneracies. The model 
furthermore provides a tool for the description of rotational bands with anomalously large gaps between the 
ground state and its first excitation. 

1 Relationship between symmetry and degeneracy: Introductory remarks. 

The relationship between the rotational symmetry of the central potentials and the degeneracy in their spectra 
is a complex one. One expects that due to rotational symmetry the spectra of the central potentials will show 
the (21 + l)-fold degeneracy with respect to the magnetic quantum number, which they really do, though not 
alone. Matters are complicated by the circumstance that various exactly solvable radial potentials (see ref. 0] for 
a compilation) have more symmetries than the rotational symmetry of the angular motion. The reason for this is 
that in the process of the reduction of the radial part of the separated Schrodinger equation down to the hyper- 
geometric differential equation by an appropriate point-canonical transformation, the energy always emerges as 
a combination, known as principal quantum number, of I, and the order n of the polynomial. In many cases the 
radial wave equation can be transformed to a power series of the Casimir operator of some Lie group different 
from SO(3) and the principal quantum number can be associated with the corresponding eigenvalues of the group 
invariant. A prominent example in that regard is the fundamental inverse distance problem, where this combination 
appears as, N = l + n + 1, the positive integer N being known as the principal quantum number. Summing up the 
(21 + 1) fold degeneracies for all I = 0, 1, 2, ..(N— 1), a larger J^ _1 (2Z + 1) = A f2 -fold level degeneracy occurs 
with respect to both I and m. The rotational (21 + l)-fold multiplicity is of course at the very root of the larger 
degeneracy. In the particular case of the Coulombic potential it is possible to associate the larger degeneracies with 
the larger symmetry group SO(4) due to the circumstance that the respective momentum space wave functions can 
be identified with the hyper-spherical harmonics on the three-dimensional spherical surface, S* 3 J2J whose isometry 
group is SO (4). Then, N — 1 = (n + I) acquires meaning of four-dimensional angular momentum. Within this 
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context, the (21 + 1) fold degeneracy is manifest at a lower level of the group reduction chain, 50(4) D SO (3) D 
SO{2). 

It is commonly accepted that in order to reveal the beyond-rotational symmetry of a central potential problem 
one has to find a suited surface of a given isometry group, Q, such that, upon an appropriate transformation, the 
radial Schrodinger wave functions can be mapped onto the bases of the irreducible representations of Q in which 
case the former would describe free geodesic motion on the Q manifold. For example, the radial Schrodinger 
equations with the esc 2 (-^), and sech 2 (^) potentials, where R and d stand for some matching length parameters 
which make the arguments dimensionless, can be transformed in their turn to the problems of the four dimensional 
rigid rotator on the three-dimensional spherical surface, 5 3 , and the hyperbolic rotator on the two-dimensional 
one sheeted hyperboloid, Hj +1 [3]. As long as the isometry groups of the S 3 - and the Hj +1 surfaces are 5*0(4), 
and 50(1, 2) respectively, the esc 2 and sech 2 potentials are said to be 50(4) and 50(1, 2) symmetric. Stated 
differently, the radial Schrodinger equations under investigation can be transformed to angular, better, rotator, 
problems on 5 3 and Hj +1 thus making their respective beyond-rotational symmetries manifest. 

The above considerations teach one two lessons. The first is that in general complete Schrodinger equations are 
not well suited as templates for studying the relationship between rotational degeneracies and symmetry because 
the rotational (21 + l)-fold multiplicity of the states may appear surmounted by the possible additional symmetries 
enjoyed by those potentials. 

The second is that the best template for studying any symmetry is a "rotator", i.e. a free geodesic motion on a 
surface whose isometry group coincides with the symmetry under investigation. 

We here are specifically interested in the relationship between rotational symmetry and degeneracy and shall 
focus attention on the rigid rotator on the two-dimensional spherical surface 5 2 whose isometry group is 50(3). 
It is described by the angular part of the Laplace operator and is known to equal (up to — 1/r 2 ) the squared orbital 
angular momentum. For a constant radius, r = a, 

V 2 = --^L 2 , r = a = const, (1) 

a 2 

the Schrodinger problem reduces to that of the rigid rotator. The quantum mechanical Hamiltonian for the free 
three dimensional rigid rotator is then the angular part of the Schrodinger equation and given by 

^L 2 Yr(6,ip) = ^l(l + l)Yr(9,v), I=na 2 , (2) 

where now a denotes the radial distance between the two bodies of the rotator, /j, stands for their reduced mass, I 
denotes their inertial moment, and Y l m (6, <p) are the L 2 eigenfunctions. These are the standard spherical harmon- 
ics, which constitute the bases of the (21 + 1) x (21 + 1) dimensional irreducible 50(3) representations, a reason 
for the observed (21 + l)-fold level degeneracy [4|. In the following we shall see that the rigid rotator problem 
provides an adequate template for studying the relationship between degeneracy and symmetry. 

To be specific, we shall perturb the rigid rotator by a cotangent function of the polar angle in which case the as- 
sociated Legendre equation can be transformed to an one-dimensional Schrodinger equation with the trigonometric 
Rosen-Morse potential. This potential is exactly solvable and one can write down exact energies and wave func- 
tions of the hindered rotator. The cotangent-perturbed rigid rotator has been mentioned independently in refs. [0, 
and Both references entertain the possibility of introducing the trigonometric Rosen-Morse potential as a non- 
central potential in the Schrodinger equation with a Coulombic potential without really executing the calculation. 
In a work (7] prior to this, the cotangent interaction has been introduced as a non-central potential in the polar part 
of the full Schrodinger problem and handled along the line of ref. 0. According to the prescription suggested in 
ref. J2j the angular equation under investigation has been subjected to a somewhat artificial transformation toward 
the hyperbolic Scarf potential with the aim to study possible 50(2, 1) symmetry properties of the latter. We here 
instead take the straight path and solve the equation for the associated Legendre functions in the presence of a 
cotangent perturbance directly. In so doing we find a result which we consider interesting to be reported on in so 
far as it sheds some surprising light on the relationship between degeneracy and symmetry. Namely, we find that 



2 



the rotational (21 + l)-fold level degeneracy of the unperturbed rotator is respected by the hindered one. 
The paper is organized as follows. In the next section we briefly review the rigid rotator problem for the sake of 
self-sufficiency of the presentation and fixing notations. Section 3 is devoted to the problem of the hindrance of the 
rigid rotator by the cotangent interaction. There, we present the exact solutions, and discuss the phenomenon of 
observing rotational (21 + l)-fold level degeneracy despite perturbance of the rotational invariance. We shall refer 
to this phenomenon as a "symmetry breaking opaqued by level degeneracies", or, shortly, "opaqued symmetry 
breaking". The paper closes with brief conclusions. 



2 The quantum rigid rotator 



The explicit expression for L 2 reads, 



1 d . n d 1 d 2 
■ sm 9— - 



sin 9 89 """89 sin 2 9 dip 2 _ 
with 9 £ [0, tt], and ip s [0, 2tt). Separating variables in eq. (fJJ amounts to 

pH(0) = l(l + l)p\ m \9), 

-m 2 <S>(p), $(y>) = e m(fi : 



1 d d m 2 
^9d9 Sm d9 ~ sh?0 
8 2 <$>(ip) 



dtp 2 

where P ; ' m (9) are the standard associated Legendre functions 0. The resulting L 2 eigenfunctions are, 



P, H (cos0)e* 



(3) 

(4) 
(5) 

(6) 



i.e. the standard spherical harmonics, Y l m (6, <p), which thereby define the basis of the (2Z+1) x (2Z+1) dimensional 
irreducible representation of the rotational group in three dimensions, 50(3). In view of the fact that eqs. © 
describe the free geodesic motion on S 2 , the standard spherical harmonics will occasionally be termed to as wave 
functions of the free geodesic motion on S 2 . 

The rigid rotator is commonly used as a tool in the description of diatomic molecules [8 |-[10|, a field that 
has recently enjoyed a significant boost through the progress in ultra-cooling techniques. The introduction of a 
potential in eq. (O has the effect to subject the rotator to an interaction, a procedure often referred to as perturbance , 
or hindrance of the rotator. Numerically solvable angular potentials like V(9) = Vq cos 2 9, and V(9) = -if (1 — 
cos 2 9) have been introduced in refs. ifTTI . and lfT2l . respectively, and subsequently employed by several authors. 

We here shall consider a rigid rotator hindered by a cotangent function of the polar angle in which case the first 
eq. (f2]i extends toward, 
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m 2 e tmip 



(7) 
(8) 



In the latter equations have admitted for the possibility that fh may not necessarily coincide in value with the 
magnetic quantum number m in eq. (0. Our case is that eq. is exactly solvable and although the cotangent 
function ceases to commute with L 2 , the interaction retains unaltered the level multiplicity characterizing the free 
rotator. Nonetheless, the system is not rotationally invariant, the signature for this being the difference between 
F(9) from eq. ® and P 1 ™ 1 (9) from eq. ©. 
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3 The cotangent-hindered rigid rotator 



In order to solve eq. Q we first change variable to 



Next we observe that substitution of eq. (O into eq. O and subsequent multiplication of both sides by Vsin^ 
amounts to the following one-dimensional Schrodinger equation (in dimensionless units) with the trigonometric 
Rosen-Morse potential, 1 / rm(^)' 

^ 1 -V RM (0)U(9) + (e+\)U(6) = 0, 

711 ( 111 -\- 1 ) ! 

v RM( d ) = \ n 2 9 -2b cote, m=|m|-->0. (10) 

The trigonometric Rosen-Morse potential is among the well known exactly solvable potentials managed by the 
super symmetric quantum mechanics (SUSY-QM) and listed in [ 1 1. We have studied eq. ( fTOb earlier in ref. lf]~3l 
and solved it independently in terms of the real non-classical Romanovski polynomials lfT4l (reviewed in ref. |[T5l ) 
instead of the Jacobi polynomials of purely imaginary arguments and complex conjugate parameters preferred by 
QJ. The spectrum of this potential problem is well known and given by lfl3l . 



e t - - \ 1 +t(t + l), t- d =n+\m\ =0,1,2,..., |m| = 0, 1, ...t. (11) 
t(t + 1J + 4 

Observing that i can take any non-negative integer value, and that |m| = 0, 1, ...t, we conclude that the < and |m| 
quantum numbers obey same branching rules as angular momentum and magnetic quantum number. In terms of 
these quantum numbers, the expressions for the wave functions are found as, 



-(t-A) 

F{9)^Ff l \e) 1 U{e)^uf A (9), Uf n \0) = N tl m\ sm t+ ^ 9e ^ Rn r ' (cot0), (12) 



where i?"'"(cot 6) denote the Romanovski polynomials, and we have attached the indices t and fh to F(0) and 
U(6). The Romanovski polynomials (reviewed in ref. 03]) satisfy the following differential hyper- geometric 
equation: 

(1 + x2) ^~ + 2 ( 2 + ^) - " (2/? + n ~ 1)R ™' = °- (13) 

They are obtained from the following weight function, 

u a ^(x) = (1 +x 2 )' 3 ~ 1 exp(-acot~ 1 a;), (14) 
by means of the the Rodrigues formula: 

1 H™ 

The parameters a and f3 for the case under investigation are calculated as, 

P = -(n+\m\ + h + l. (16) 



n+|m| + | 2 
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Figure 1: The splitting, AE = (E t — E t -i), for 6=1, between neighboring rotational levels belonging to the 
ordinary rigid rotator (circles connected by a dashed line) and the cotangent-perturbed rotator (squares connected 
by a solid line). Noticeable, the anomalous splitting between the ground state and its first excitation. For the higher 
lying excitations the perturbed level splittings rapidly approach the unperturbed ones, a behavior visualized by 
practically coincident squares and circles. 



The labeling of the wave functions [/„ (9) corresponds to, 



Ur\6) =iV tH e-^(l + cot 2 



K£t-\fh\( cot 



with a and (3 from eq. ( TT6b . 
The functions 



Xr(0,<p)=FFK0)e 



imtp 



i in 



(17) 



(18) 



\/sm9 

determine the solutions of the cotangent-perturbed rigid rotator on S 2 . 
In summary: 

• The branching rules between the non-negative integers t and \m\ quantum labels introduced in eq. (fTTT i are 
indistinguishable from those of orbital angular momentum and magnetic quantum number. In this sense the 
energy spectrum, e t , in eq. ( fTTb formally shows same (2t + 1) fold level multiplicity as the unperturbed 
rotator. The splitting between two subsequent levels reads, 



Et — Et-i — 2t 



2tb 2 



2t + 



(19) 



Equation $1% shows that for b ~ 1 the perturbation affects mainly the splitting between the ground state and 
its first excitation and becomes rapidly negligible afterward. In this fashion the rigid rotator hindered by a 
weak cotangent interaction describes a rotational band with an anomalously large gap between the ground 
state and the first excited level according to E\ — Eq = 2 + ^-b 2 (see Fig. 1). 
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Figure 2: The dependence of the damping factor (i.e. the exponential) in eq. (fT2l on i. The flattening with the 
increase of t (from bottom to top) is well pronounced. 

• Though the X" l (9, ip) functions in eq. ( fl~8b do not transform according to a irreducible representation of L 2 
in eq. (0, for large t values they asymptotically approach the spherical harmonics according to: 

sin* 0PT 3 ' 2 (cot0) ^ sin* 9 R^'}^' (cot 9) ~ P t |m| (cos 0). (20) 
This is due to the rapid flattening of the exponential (damping) factor with the increase of t (Fig. 2). 
The last equation in the chain has already be found in ref. [7 1 from a different perspective. 

It is quite insightful to take a closer look on the structure of the P t ' m '(6*) = t/J m ' (9) / V sin 9 functions in 
eq. (fT2l . Below we decompose some of the lowest functions into associated Legendre functions. In omitting the 
normalization factor for simplicity, we find: 

n = 1 : fh = : Ff(9) = e -2 ^ sin 6nf ~ 3 (cot 9) = er^ (-P?(cos6) + yP^cosfl)^ , (21) 

n = l: m = l: F^(9) = smORj '"^(cot6») = e^P^cosfl), (22) 

n = 2: m = 0: F 2 °(6») = e -3 ** sin 2 flPj^ 1 (cot 6») = e~™ (ip 2 °(cos6») + ^(P^c 

+ jjp 2 (cos#))), (23) 

n = 1 : m = 1 : ^(0) = e _a r s in 2 0pf ~* (cot (9) = e -3 ^ (P^cosfl) - ^P 2 2 (cos(9)) , (24) 

9,„n 268 9 M 



cos ( 
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n = 0: m = 2: P|(0) = e - ^ sin 2 9RJ '~*(cot8) = e~^P 2 2 (cos (9). (25) 



The latter equations show that the wave functions f\ m ' (0) of the cotangent-hindered rotator are linear combinations 
of damped associated Legendre functions, P ; ' m ' (cos 9), with I = t and m € [m, i]. The transformation of the wave 
functions of the perturbed problem to the wave functions of the free rotator is clearly non-unitary and makes the 
breaking of the rotational symmetry at the level of the representation functions manifest. The above considerations 
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provide a broader view on level degeneracy compared to the algebraic methods |[T6l . which conjecture that the 
symmetry of the Hamiltonian is same as the symmetry giving rise to the observed degeneracies in the spectrum. 
Specifically, the algebraic Hamiltonian of any type of a (21 + 1)- fold level multiplicity is supposed to be described 
by means of a Hamiltonian that is a power series of L 2 . Specifically the algebraic Hamiltonian, W, for the level 
multiplicity in eq. (TT2l i would be (in dimensionless units): 



L 2 + ±' 



(26) 



effectively equivalent to 



H 



l 2 (l + lf 



4(i(Z+l) 



(27) 



in which case the standard spherical harmonics would remain eigenfunctions of the perturbed rotator, 



1(1 + 1)- 



1(1 + 1) + 



Yi m (0^)- 



(28) 



The difference between the algebraic method in eqs. (|26>-(1281. on the one side, and the Hamiltonian method in 
eqs. ( [Tol l. (fT2l) . and ( fT8l will become detectable at the level of transition probabilities. 



4 Conclusions 

In this work we considered the option of having a symmetry breaking at the level of the representation functions 
while preserving the degeneracy patterns in the spectra of the original symmetry. It this fashion, we argued that a 
symmetry breaking opaqued by level degeneracy can take place. On the example of a cotangent-perturbed rigid 
rotator we showed that the (21+ 1 )-fold level degeneracy is a softer criterion for rotational symmetry realization, the 
rigorous one being the property of the wave functions to transform according to irreducible representations of L 2 
in (01. In eqs. (F2TI) -(f25]) we noticed that the solutions, F(6), of eq. (0 can be represented as linear combinations 
of (exponentially) damped associated Legendre functions. In terms of the latter, damped spherical harmonics, 
Y™(0, (p), can be defined according to, 

e-^L 2 e^Y t m (6,<p) = 1(1 + l)Y t m (6,tp), 

Y t m (9,(p) = e~^P t m (9)e mv (29) 

In Figs. (30 we compare the regular and damped spherical harmonics for illustrative purposes. 

The similarity transformation of L 2 in the last equation is a dilation transformation which shows that the free 
and the hindered rotator solutions are related by a non-unitary transformation. The equations (l2H -(l25t in combi- 
nation with d29l allow to give the cotangent-hindered rotator also the reading of a realization of an so(3) algebra 
unitarily inequivalent to the standard rotational. Specifically for the case of a weak perturbance, we revealed possi- 
bility of observing in diatomic molecules rotational bands characterized by an anomalously large gap between the 
ground state of the rigid rotator and its first excitation, while practically preserving the subsequent level splittings. 
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Figure 3: The regular spherical harmonic |Yj°(0, f)\ (left) versus the damped, \Yi(9, <p)\, (right). 
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Figure 4: The regular spherical harmonic |Kj (#, <p)\ (left) versus the damped, |i^ 2 (^i f)V (right). 
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